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One-electron self-energy in the t-J model was computed using a recently developed large- N 
method based on the path integral representation for Hubbard operators. One of the main fea- 
tures of the self-energy is its strong asymmetry with respect to the Fermi level, showing the spectra 
mostly concentrated at high negative energy. This asymmetry is responsible for the existence of 
incoherent structures at high negative energy in the spectral functions. It is shown that dynamical 
non-double-occupancy excitations are relevant for the behavior of the self-energy. It is difficult to 
understand the asymmetry shown by the self-energy from weak coupling treatments. We compare 
our results with others in recent literature. Finally, the possible relevance of our results for the 
recent high energy features observed in photoemission experiments is discussed. 
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It is commonly accepted that high-T c cuprates are 
strongly correlated systems. In these materials angle 
resolved photoemission spectroscopy (ARPES) experi- 
ments show interesting and still unexplained features: a) 
the low energy kink^ at energy ~ 40 — 70meV and; b) 
the high energy anomalies^^^ at ~ 0.5 — leV, known 
as waterfall. For a theoretical description of these ex- 
periments it is necessary to calculate self-energy correc- 
tions on electronic correlated models as t — J or Hub- 
bard. In spite of the progress by means of numerical 7 - 
and analytical^ methods, the problem remains of huge 
interest. Recently we have proposed a large- N approach 9 
for the t — J model which is based on the path integral 
representation for Hubbard operators (called PIH in what 
follows) . In the PIH method the spin index a is extended 
to a new index p running from 1 to A" and the pertur- 
bation is performed in powers of the small parameter 
1/N. It was shown that in leading order of 1/N (0(1)), 
which is equivalent to mean field level, PIH results agree 
with the slave boson& calculation. The large- N expansion 
provides a controllable way for selecting and truncating 
Feynman diagrams. However, the results are more rep- 
resentative for the physical case N = 2 when terms in 
powers of 1/N can be collected. In this context, PIH can 
be implemented beyond mean field allowing the calcula- 
tion of self-energy corrections and spectral functions. 10 ^ 1 
The obtained spectral functions were compared with ex- 
act diagonalization results finding good agreementj 1 ^ 11 

In this paper we discuss the role of dynamical fluctua- 
tions of the non-double-occupancy constraint on the self- 
energy results. In addition, differences between present 
self-energy and that obtained from calculations based 
on weak coupling approaches like, for instance, random 
phase approximation 12 (RPA) are discussed. We com- 
pare also our results with those obtained by other cal- 
culations, and discuss the possible relevance of present 
results for the recent high energy features observed in 
ARPES experiments in cuprates. 

In Refs. 

sua 

it was discussed that PIH approach weak- 



ens collective spin fluctuations over charge fluctuations. 
Although for finite doping away from half filling the rele- 
vance of magnetism is a matter of debate, for preventing 
possible objections about the influence of magnetic con- 
tributions, we calculate for the high doping value 5 = 0.3. 
This high doping corresponds to highly overdoped regime 
of cuprates where magnetic fluctuations are found to be 
very weakJ^ In addition, and for simplicity, we present 
results for J = 0. For high doping, PIH does not show 
strong dependence with J, being representative the re- 
sults for J = (see Ref.[14j] for discussion). On the other 
hand, no strong J dependence is expected for high doping 
values. 

Collecting all 0(1/ A) contributions, the full self- 
energy S(k, ui) (real and im aginary parts) in the square 
lattice is described in Refs. jld l llj |. Herein, for conve- 
nience, we reproduce only the corresponding results for 
scattering rate: 

+ ImS AA (k, W ) (1) 

where 

Im Eiyj(k, u>) = Im[D RR (q, u - Sk- q )} 

6 q 

X [n F {-£k-q) + n B {u - £k-g)] > 
JmSflA(k,u) = — - H Im[D\ij(<i, u ~ £fc-g)] 

1 * s 

q 

x [n F (-e k - q ) + n B {oJ - e k - q )\ , 
ImY,\\(k,ui) = — Im[D\\ (q, ui ~ Sk- q )] 

1 * s 

q 

X [n F (-£k- q ) +n B {u - e k - q )] . (2) 

In eq.|2]), £1 = (efc_ 9 +cj + 2/i)/2 and e k — —tS(cosk x + 
cosfcj,) — /i is the mean field electronic band. N s is the 
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FIG. 1: (Color online) Spectral functions At (solid lines) and 
Arr (dashed lines) from F (panel (a)) to ki? (panel (f)) in 
the nodal direction for 5 = 0.3. Energies are in units of t. 
The vertical line at lu = marks the Fermi level. 



number of sites, np (^b) is the Fermi (Bose) factor and 
^ the chemical potential. Djjn(q 7 lu) is the charge-charge 
correlation function which contains collective charge ex- 
citations. D\\ and Dr\ correspond to the pure non- 
double-occupancy sector and the mixing between non- 
double-occupancy and charge sectors respectively. 

Starting from Y,t and X.r.r the spectral functions 
Ar(k, w) and ArrQs., lu) are respectively calculated. In 
the calculation of the spectral functions, the real part of S 
was numerically computed by using the Kramers-Kronig 
relation from eq.(l). Spectral functions Ar(k, cj) (solid 
lines) are presented in Fig. 1 from T (k = (0,0)) (panel 
(a)) to the Fermi vector kp (panel (f)) in the T—(tt, tt) di- 
rection (nodal direction) of the Brillouin zone (BZ). The 
sharp peak near u) = is the quasiparticle (QP) peak 
and defines the QP Zhang- Rice 15 (ZR) band of the t — J 
model whose bandwidth is reduced from that of the mean 
field band due to the self-energy renormalizations. In 
addition Ay(k, to) shows incoherent spectrum (IS) which 
is mainly localized at high negative energy (lu ~ — At). 
Almost no IS is observed for lu > 0. Similarly to Lanc- 
zos diagonalization results ) 10 ! 16 ' 17 while the QP disperses 
through the Fermi surface, the IS moves in opposite di- 
rection. Present results are for temperature T = OK, and 
a finite value for T does not change the main conclusion. 
As discussed below, self-energy contributions (Fig. 3 ), re- 
sponsible for the IS, lie on an energy scale of the order 
of several t. Therefore, no significant changes occur for 
realistic values of temperature. 

Fig. 2a shows, in the main directions of the BZ, the 
energy dispersion of the QP peak (solid circles) and the 
IS (open circles) shown by At- In panel (b) we repro- 
duce the t — J model results from Fig.lc of Ref.[l8[ ob- 
tained, also for doping 5 = 0.3, using Gutzwiller pro- 
jection variational Monte Carlo (VMC). Although both 




FIG. 2: (Color online) (a) Energy dispersion of the QP peak 
(solid circles) and the IS (open circles) shown by At for 5 = 
0.3. (b) Results reproduced after Ref.pl| for a qualitative 
comparison with panel (a). Dotted line in both panels marks 
the Fermi level. 



methods (PIH and VMC) are different, results are sim- 
ilar. Both panels show a QP (or ZR) band and IS at 
high negative energy. None of the obtained results show 
significant IS at positive energy. Although the IS is more 
dispersing in PIH than in VMC, its energy position is 
of the same order of magnitude in both methods. The 
QP band is somewhat flatter in PIH, if t = 0.4eV.-^ For 
instance, the QP at T is at lu ~ —OAeV for VMC and at 
lu ~ —0.2eV for PIH. (In Ref.[10( it was discussed that 
the QP bandwidth predicted by PIH is reduced from that 
obtained by Lanczos). The size of the circles scales lin- 
early with the spectral weight (SW). The SW on the BZ 
is more homogeneous in panel (a) than in panel (b). PIH 
predicts, at kp, a QP weight Z ~ 0.5 indicating that, 
even for S = 0.30, ~ 50% of the SW is concentrated in 
the IS. 

The scattering rate — /mS-rfk, to) at k = kp in the 
r — (tt, tt) direction is presented in Fig. 3. JmEy (solid 
line) is very asymmetric with respect to lu = showing 
most of the SW at lu < 0. This asymmetric behavior 
is the cause of the shape of At in Fig.l. This strong 
asymmetric distribution should be interpreted as a con- 
sequence of the difference between addition and removal 
of a single electron in a correlated system. Recently^ us- 
ing Lanczos diagonalization in the t — J model, a similar 
asymmetric behavior was also discussed. It is important 
to notice that our scattering rate does not show any low 
energy scale, thus it can not explain the low energy kink. 
If the low energy kink is due to magnetic excitations, 20 
or other electronic effects, they are obviously not include 
in our approach. However it is possible that the kink is 
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FIG. 3: (Color online) Scattering rate for Er (solid line) and 
E-RR (dashed line) versus uj for S — 0.3 at k = ki?. Inset, 
scattering rate predicted by RPA on the Hubbard model for 
8 = 0.3 at kF (see text). 

due to phononsi^i 

For understanding which are the leading contributions 
responsible for the asymmetric behavior of St, we have 
calculated —ImT^RR- In contrast to /toS^, Im^RR 
(dashed line in Fig. 3) is very symmetric with respect to 
uj = and App (dashed lines in Fig.l) shows the IS dis- 
tributed almost equally at positive and negative uj. Evi- 
dently S^a and S^a (eq.(P) and eq.(2)) contribute signif- 
icantly, showing the relevance of non-double-occupancy 
excitations on the redistribution of the SW for leading 
to the final form shown by At- In spite of the doping 
studied here corresponds to the highly overdoped regime 
of cuprates and it could expected that a weak coupling 
approach, like RPA, be reliable, this is an open and con- 
troversial issued For instance, and important for present 
discussion, ARPES experiments show, for highly over- 
doped samples, high energy anomalies with similar char- 
acteristics to those in underdoped samples ,SiMiSiS It is 
worth to mention that £_r.r (first line in eq.(2)) has a 
somewhat similar meaning to the self-energy when only 
charge fluctuation are considered in RPA for the Hub- 
bard model. This later can be written as- 3 

ImE RPA (k,oj) = --L^([//2) 2 /m[ Xc (q,w-£fc- g )] 

S q 

x [n B (oJ -&-,) +nj?(-£fe-g)] 

(3) 

where Xc is the RPA charge susceptibility^ JJ is the on- 
site Coulomb repulsion, and = —2 i(cos k x +cos k y ) — ^ 
is the bare tight-binding band on the square lattice. Since 
in both, ImT, pp and IrnSppA, charge excitations are in- 
volved in the electron renormalizations, it is instructive 
to compare results from both sides. Inset in Fig. 3 shows 
—ImYjRPA at the nodal kp for U = 4 and 8 = 0.3. 



Interestingly, ImY^ppA has a similar shape to that for 
ImSjjii, i.e., —ImYiRPA is very symmetric with respect 
to uj = 0, leading also (not shown) to IS homogeneously 
distributed at both sides of the Fermi level as for the 
case of App (Fig.l). We have found this behavior for 
Erpa very stable against different conditions for U and 
hole doping away from half-filling. T,pp and UppA show 
similarities because they can be simply interpreted in 
terms of fermions interacting with charge fluctuations. In 
our opinion this self-energy behavior, symmetrically dis- 
tributed around uj = 0, can be expected in weak coupling. 
However, results are different when the full self-energy 
(St) is considered. St is obtained in strong coupling 
and contains fluctuations above mean field level which 
are very difficult to be obtained perturbatively from usual 
fermions. The strong coupling calculation suggests that 
electrons interact with charge fluctuations and with ex- 
citations which represent non-double-occupancy effects. 
These excitations, expressed in our approach by Dp\ and 
D\\, are responsible of the concentration of the incoher- 
ent spectral weight at negative energy. In addition, they 
are dynamical (q and ui dependent) and, they can not be 
simply considered as a static enforcement of non-double- 
occupancy constraint as in mean field approximation. At 
this point we wish to emphasize about the important 
role of the non-double-occupancy constraint even for the 
high doping studied here. In addition to the results in 
Ref.[I| for 8 = 0.30, Lanczos results^ for 8 = 0.25 show 
also large IS at negative uj. Since this behavior can be 
understood if the scattering rate is asymmetric with re- 
spect to uj = 0, we think that these results indicate that 
the overdoped t — J model shows strong coupling fea- 
tures. With increasing doping, our results will be closer 
to those obtained using RPA. For 8 > 0.7, St becomes 
almost symmetric and, in this case, T,pp approaches St 
showing that S^a and Saa have little influence. 

Next, we discuss the possible relevance of present re- 
sults for the high energy anomalies observed in ARPES 
experiments in cuprates. Momentum distribution curves 
(MDC) analysis of the experiment suggests the occur- 
rence of one band which is strongly renormalized near the 
Fermi surface. Away from the Fermi surface this band 
develops an abrupt change reappearing at high energy 
(~ —leV), given the impression of a waterfall. This re- 
sult supports a description in terms of renormalizations 
of the LDA 3 -^ or an uncorrelated bandit In contrast, 
energy distribution curves (EDC) analysis shows the si- 
multaneous presence of both, low energy and high energy 
excitations^ These results suggest the occurrence of a 
low energy band, associated with the ZR band of the t — J 
model, and IS at high binding energy 18 i 19 ' 25 Therefore, 
our results are in closer agreement with the interpreta- 
tion obtained from the EDC analysis. At this point it is 
important to remark that it was recently discussed 2 ^ that 
the waterfall dispersion is not an intrinsic feature but re- 
sults from the suppression of the photoemission intensity 
near r due to matrix elements effects (see also Ref. [6] for 
discussion about the momentum and energy distribution 
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curves dichotomy). Concluding, we propose that the high 
energy features can be described in the framework of the 
t — J model which shows the existence of a low energy 
ZR band and IS at high negative energy. Additionally, 
we have shown that dynamical non-double-occupancy ex- 
citations are relevant for transferring most of the SW to 
high negative energy, leading to a well pronounced IS at 
lj < as the observed by the experiment. In addition 
to the requirement that the IS should be mainly concen- 
trated at negative w, its SW should be large enough to 
be observed. As discussed above, this condition is also 
satisfied by the t — J model. 

In summary, we have discussed that dynamical non- 
doublc-occupancy effects, which are only obtained be- 



yond mean field level, are responsible for a strong asym- 
metry of the self-energy with respect to d = 0. This leads 
to spectral functions where large IS is present at high neg- 
ative energy with nearly no signals of IS at positive u>. 
It was also discussed that this picture is very improba- 
ble to be obtained from methods which treat the elec- 
tronic correlations in weak coupling. Our results show 
similarities with the recent high energy features observed 
by ARPES experiments in cuprates giving an additional 
support to the point of view that electronic correlations 
push cuprates to the strong coupling regime. 
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